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On the Mean Duration of Random Walks 1 

Wolfgang Wasow 

The mean duration of a discrete random walk in a bounded domain is studied for 
general transition probability functions. The mean duration satisfies an integral equation 
which, for small mean step length, can be approximated by an elliptic differential equation. 
This leads to explicit inequalities for the mean duration in an important special case. 



I. Introduction 

This paper is concerned with random walks in a 
bounded open domain G in n -dimensional Euclidean 
space. The walk ends when the moving point leaves 
G for the first time. The nature of the random 
walk is characterized by the transition distribution 
function F (Q, P). Here P and Q are points of the 
space, and F (Q, P) is the probability that a point 
at P will, after one step, be in the rectangular sector 
of the space in which all coordinates are less than 
those of Q. 

By duration of a random walk we mean the total 
number of steps of the walk, including the one on 
which the point leaves 67. The subject of this paper 
is the mean duration; that is, the expected value 
W (P) of the duration as a function of the starting 
point P. 

It will be shown that — under very general condi- 
tions — W (P) is the unique solution of the integral 
equation (9). In section III assumptions are 
introduced, which, in a sense to be specified below, 
make the occurrence of all but very small single 
steps very unlikely. It will be shown that under 
those assumptions W (P) differs very little from the 
solution of a certain elliptic differential equation. 
This differential equation involves only the first and 
second moments of F (P,Q). The method used in 
these sections is closely related to that of Petrovskv 
[4]. 2 

For a particularly important class of random 
walks this differential equation reduces to a special 
case of Poisson's equation. For such random walks 
it is possible to obtain numerical inequalities for the 
mean duration. These estimates confirm the plausi- 
ble conjecture that the mean duration does not grow 
very fast with the number of dimensions, a fact 
that is of interest in connection with the attempts 
to use stochastic sampling procedures for the 
numerical solution of partial differential equations. 

If the moving point performs a continuous 
Brownian motion rather than a discrete random 
walk, the mean duration satisfies exactly a Poisson 
equation. This is an easy consequence of the fact 
that the fundamental solution of the differential 
equation for diffusion in the presence of absorbing 
boundaries can be interpreted as the probability 
density that a Brownian particle starting from a 
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given point P at time be at a point Pi at time t 
without having met the boundary. (See, e. g. 
Smoluchovski, Drei Vorlesungen Uber Diffusion, 
etc., Physikalische Zeitschrift 17, p. 568 to 570 
(1916), or P. Levi, Les processus stochastiques, 
p. 272 to 273 (Paris, 1948)). It has been pointed 
out to the author by M. Kac that this fact could be 
used as a starting point for an alternative treatment 
of our problem based on the theory of probability 
measures in function space. 

In sections I to III, where the w-dimensional case 
differs from the two-dimensional one only by the 
need for more cumbersome notation, we shall give 
the arguments for two dimensions only. 

II. Preparations; Lemmas on Integral 
Equations 

The transition distribution function F (Q, P) is 
assumed to be bounded and Borel measurable for 
all P and Q and a distribution function with respect 
to Q. In order to guarantee that the moving point 
leaves the domain G with probability one we intro- 
duce a hypothesis which will be called assumption 
(A). 

Assumption (A). There exist two positive num- 
bers r and s, independent of P, such that the proba- 
bility is at least s that the abscissa of the moving 
point be increased by r or more in one step. 

We refer to [4, p. 431] for the simple proof that 
the probability of leaving G is indeed one, if assump- 
tion (A) is satisfied. 

Denote by F m (Q, P), m=l, 2, . . . , the iterated 
functions defined by 



F 1 (Q,P)=F(Q,P) 



FUQ,P> 



■//■ 



.,((2, R)dF(R, P). 



Here, and in the sequel, the Stieltjes differentials 
are formed with respect to the first argument, and 
an integral without explicitly indicated domain of 
integration is to be extended over the whole space E. 
F m (Q, P) is the probability of a point being in the 
domain x<C%q, V^Vq after exactly m steps, starting 
from P and having all of its intermediate positions 
in G. Here (x Q , y Q ) are the coordinates of Q. In 
consequence of assumption (A), and the bounded- 
ness of G, there exists a positive integer m and a 
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positive number c, less than one, such that 
( dF m (Q,P)<c, for all P in 6. 



(1) 



We shall need the following lemma concerning 
integral equations formed with distribution func- 
tions. For a proof, we refer again to [4]. 

Lemma 1: If assumption (A) is satisfied, there 
exists exactly one bounded function u(P) which 
solves the integral equation problem 



u(Py- 



fu(Q)dF{Q, P) in G 
i(P) outside G, 



(2) 



where \p(P) is a prescribed bounded and Borel meas- 
urable function. This solution can be found by 
means of successive iterations starting from any 
function that is equal to \[/(P) outside G, and bounded 
and Borel measurable in G. 

Corollary: The solution u(P) of (2) satisfies the 
inequalities 



g.Lb.^(P)<u(P)<l.u.b.^(6). 



(3) 



Proof oj the corollary: Take as initial approxima- 
tion to u (P) the function u (P), which is equal to 
g.l.b.iKP) in 0. Then 

um=fuo(Q)dF(Q,P)>g.\.b.f, 

and, by induction, u n >g.\.b.\p. A passage to the 
limit proves the first one of the inequalities (3). 
The proof of the second inequality is analogous. 
Lemma 2: Let </>(P) be nonnegative and Borel 
measurable in G, and assume that assumption (A) 
is satisfied. Then the integral equation 



u(P)= f a u(Q) dF(Q,P) + 4>(P) 



(4) 



has a unique bounded solution in G. If <t> n {P) is a 
nondecreasing sequence of nonnegative functions 
converging to <t>(P) , and u n (P) is defined by 



u o =0,u n aP) = j G Un(Q)dF (Q,P) + $ n (P), (5) 



then 



Z7(P) = lim u n (P) 



(6) 



is the solution of (4). The solution u(P) is non- 
negative. 

Prooj: The sequence u n (P) is proved ;to be non- 
decreasing by means of mathematical induction, 
using the relations 

*i (P)-tio(P)= f G ^(Q)dF(Q,P)>0 3 
and 



u n+1 (P)-u n (P)=j G hn (QD-UnMQ)] dF (Q, P)+ 

<l> n+1 (P)-<t> n (P). 

In order to show also the 4 boundedness of the se- 
quence u n (P) we iterate (15) and find 

u n UP)= f f u n - 1 (Q)dF(Q,R)dF(R,F)+ti 2) (P) 

J GJ G 

= j G n n _ 1 (Q)d Q [f/(Q, R)dF(R } P)]+ </>f> (P) 

= J G n n . 1 (Q)dF 2 (Q,P) + ^ (P), 

where 

A? (P)=j G <l>n-i(Q)dF{Q,P)+<l> n (P)<2 l.u.b.*(F). 

Repeating this procedure m times — assuming n>m 
— the integral equation 



U n 



ln + l(P) = j Q U n - m (Q)dF m (Q 1 P)+<t>i n+l) (P) 

is obtained, where 

flJ»+»(P)<(m+l) 1. u.b. </>(P). 

If L n denotes the 1. u. b. of u n (P) in G we have 
therefore, using (1), 

L n+ i<cL n - m + (m + l) 1. u. b. </>(P) 



hence 



<cL n + (m+l) l.u. b. </>(P), 



L n <(l-c)- l m 1. u. b. 0(P). 



The sequence u n (P) , being bounded and nondecreas- 
ing has therefore a limit, and this limit satisfies the 
integral equation (4), since a passage to the limit 
under the integral sign in (5) is legitimate. 

The uniqueness of the solution follows by a stand- 
ard argument from the uniqueness of the solution 
of the corresponding homogeneous problem. The 
uniqueness of the latter is assured by lemma 1 . 

Corollary: Equation (4) has a unique solution, 
even if the assumption that </>(P) be nonnegative is 
omitted. 

Prooj: Let ^ be a negative lower bound of <£(P). 
Lemma (2) carTbe applied to the integral equations 

u>(P)=f G Ui(Q)dF(Q,F)-* 

u 2 (P)=j G U 2 (Q)dF(Q,P) + <t>(P)-£ 
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u(P)=u 2 (P)—u } (P) is then the solution of (4). 
Lemma 3: If </> and <j> ar c constants such that 

±<4>(P)<*volG, 

and if W(P) is the solution of 

W(P)=j G W(Q)dF(Q,P)+l, 

then the solution of (4) satisfies the inequalities 

±W(P)<u{P)<l>W(P). (7) 

Proof: The function u(P) =<£ W(P) is the solution of 

u(P)= f u(Q)dF(Q,P) + $. 

JG 

Subtracting (4) from this we see that u{P)—u{P) 
is the solution of a similar integral equation with 
0— <£(P) as nonhomogeneous term. From lemma 
2 it follows that u(P)—u(P) is nonnegative. The 
other inequality of (7) is similarly proved. 

Ill, The Integral Equation for the Mean 
Duration 

Let p n (P) be the probability that the duration of 
a random walk starting at P be exactly n. Because 
of the remark following the assumption (A) , we have 



S?.(f)=i, 



(8) 



if we assume — as we shall — that assumption (A) is 
satisfied. 

For a point outside G 

Pn(P) = d 0n (P not in 67) 

where 5 0n is Kronecker's symbol. The expected 
duration of a random walk is — by definition — 



Let 



W{py 



■■Hn Vn {P) 

71 = 



if this series converges. Otherwise, the expected 
duration is infinite. We shall prove 

Theorem 1. If assumption (A) is satisfied, W(P) 
is finite and is the unique solution of the integral 
equation 



Proof: 
tion 



W(P) = j G W(Q)dF(Q,P) + l. (9) 

The probability p n {P) satisfies the rela- 

p n+ i{P) = jpn(Q)dF(Q,P),n>0. (10) 



W n (P)=Jlv Pv (P), 
then we obtain from (10) the relation 

W n+l (P) = fw n (Q)dF(Q,P) 



71+1 
p=0 



Reference to lemma 2 completes the proof. 

If F (Q, P) is a discrete distribution with a finite 
number of discontinuities, and if the position of these 
discontinuities is independent of P, then the integral 
equation (9) reduces to a difference equation. The 
simplest cases can be solved explicitly; see, e. g. 
[3, p. 286.] 

IV. The Asymptotic Differential Equation 

In many applications the steps of the random walk 
are likely to be small. In order to describe more 
precisely the sense in which this statement is to be 
understood we assume, following Petrovsky, that the 
transition function depends on a small parameter 
ix in such a way that the following assumptions are 
satisfied. 

Assumptions (B). Denote by a t (P, /z), b ik (P, /*), 
(i, k=l, 2) the first and second moments of F (Q, P) 
about the point P. Then 

a,i(P,fx) = ai(P)tx+o(ix) 

b ik (P,ix) = P ik {P)ix + o{ix). 

The functions o (jj.) may depend on P, but they are 
to have the indicated order uniformly in G. The 
a t (P) and p ik (P) are assumed to be in class C i2) in 
G+O. 

Assumption (L). Let K r (P) denote the circle 
P0<r, with center at P. Then 



! 



P$dF{Q,P)=o{p), 



E-Kr(P) 

here o (fx) depends on P and r. The relation is to 
hold uniformly for P in the closure of G, for any 
fixed r>0. 

In the applications ll may measure the smallness 
of the mesh, if the walk takes place in a lattice, the 
smallness of the standard deviation, if the transition 
distribution is normal, the smallness of the time 
intervals between observations, if the walk is actually 
a motion in time, etc. 

The definition of the fi ik is readily seen to imply 
that jSii 022— 0i2>O. We require, beyond this, 

Assumption (E). In G and on its boundary, 

0ii022-/3? 2 >O. 

This assumption will enable us to make use of the 
existence theorems in the theory of elliptic differential 
equations. 
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Finally, we have to require a certain degree of 
smoothness of the boundary C of G. In particular, 
our results remain valid for polygonal domains. 

Assumption (S). C has a continuously turning 
tangent. 

This restriction can be weakened. The necessary 
arguments can be found in [4, p. 438]. 

We shall show that the function w = nW remains 
finite, as /x— >0 and converges to the solution of the 
differential problem 



V=0 onC 



(11) 



where, 



&[V]=lfi u V a +fi u V, t +±PnV n +a l V.+a t V,. 

(12) 

Because of assumption (E) the differential problem 
(11) is elliptic and possesses in G a unique solution. 
The basic idea of the proof is to show that V 
satisfies an integral equation little different from the 
one for W. In order to avoid extraneous difficulties 
near the boundary C it is convenient to replace 
V(P) by a slightly different function Vs(P) defined 
in a larger domain G n . This can be done by con- 
structing a twice continuously differentiable mapping 

x'=f(x, y;S), y'=g(%, y;S), (x, y in 67) 

which is continuous in <5, for <5>0, together with its 
first and second derivatives with respect to x and y, 
and has the following properties: (a) It reduces to 
the identity for 8— >0. (b) It is, for all 5, the iden- 
tity transformation in a subdomain G' of G that 
tends to G as $—>0. (c) It maps G into a domain G" 
containing G in its interior for <5>0. 

For the explicit construction of such a mapping 
with the help of assumption (S) we refer to [4] . 

If we define V 8 (P) in Q by 

V b (x',y') = V{x iy ),{x,ymG) 

then this new function is defined and twice con- 
tinuously differentiable in G". It tends to V(P) y 
uniformly in G together with its first and second 
derivatives. We extend the definition into the whole 
plane E by setting 

V s (P) = in E-G". 

In order to estimate the integral 



J=jv i (Q)dF(Q,F) 



for P in G we denote by K(P,8) a circle about P 
with radius r(8) depending on 8 but not on P, such 
that K(P,8) is in G" whenever P is in G. We define 



J,= ( V i dF 9 J 2 =[ V 8 dF, 

JK(P,8) JE-K(P,8) 

then J=J x -\-J 2 . 

Now if A is an upper bound of | Va|, for all 8 with 

0<<5<5 , 

|</ 2 |<-4f PQ 2 dF(Q,P). 



Hence, by assumption (L), 

J 2 =o (/i), for fixed 8. 



(13) 



For the calculation of J x we use Taylor's formula, 
for PmG,Q in K(P, 5), setting 

V 8 (Q)=V 8 (P)+M[V s (P)] + Q[V d (P)]+R[V 8 ], (14) 

where 

M[V] = (£-x) V x (P) + (v-y) V V (P) 

Q[V]=\{^-x)W xx {P) + 2^-x){y ] -y)V xv {P) + 



and 



i. e., 



with 



( v - y yv m (P)} 



R[V] = Q[V(P)]-Q[V(P)],PmK(P, 5) 



iW]=PQ 2 P (P,Q, 8) 



(15) 



lim p(P, Q, 5) = uniformly for P in G, Q in K(P, $). 

(16) 

In these formulas x, y and £, 77 are the coordinates 
of P and Q, respectively. We now integrate (14) 
termwise with respect to Q over K(P, 8) and discuss 
the terms separatelv. We have, e.g., by assump- 
tion (B), 

f (Z-x)dF= l xa x + o( i i)- f (S-x)dF, 

JK(P,S) JE-K(P,S) 



but 



\jE 



(£-x)dF 



^H8)°^' 



lE-K(P,d) 

by assumption (L), and hence 

f ^-x)dF= f ia l + o(fi), forfixed 8. (17) 

JK(P,5) 

Similar arguments apply to 
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Also, 



f dF, f (n-y)dF, etc. 

JK(P,8) JK(P,8) 



L 



K(P,8) 



R[V i ]dF= Pl (P,8)0(ji) 



(18) 



by virtue of assumption (B). Here pi(P, 5) tends to 
zero, uniformly in 67, as 5— >0, because of (1.6), 
and 0(/x)//x remains bounded as 5-»0. Formulas 

(1.3), (1.7), (1.8) and the formulas for f dF, 

J K(P,8) 

(r)—y)dF, etc., analogous to (1.7), when 



/, 



K(P,8) 

applied to (1.4), yield 

J=VW + vL[V s ] + Pi (P,8)0(ji) + o(ji). (19) 

Now let e>0 be assigned arbitrarily (but less than 
A), and choose a <5, (0<<5<<5 ), so small that the fol- 
lowing inequalities hold: 



(a) 



\L[V d ] + l\<e/2±AmG 



(20) 



This is possible because of (11) and the properties 
of Vs as described above. 



(b) | Pl (P,5)0( M )|O/24A 

Here the remark after formula (18) is used. 



(c) 



(d) 



\V s -V\<~inG 



\V i \<^mG"-G 



(21) 



(22) 



(23) 



This last inequality is the only point where the 
boundary condition in (1 1) is needed. 

After 5 has been fixed in this fashion we choose 
Mo>0 so small that the o(n) in (19) satisfies the 
inequality 



|o(ju)|</W24A in 6, for | M |<Mo- 



(24) 



From (20), (21) and (24) it follows that (19) can be 
written 



J 



V*(Q) dF(Q,P)=Vt(F)-*9<P,v), (25) 



where 



|(/(P, M )-l|<e/8AforPin^| M |<Mo. (26) 



For the comparison of Vh with pW we split Vs 
into 

where 



l 



W>dF-- 



y«'- M </,Fa (1) = 



-0 in E-G (27) 



f V? ) dF= VJ» in G, V^=V„ in E-G. (28) 

The integral equations (27) and (28) have unique 
solutions, by the corollary to lemma 2, and by 
lemma 1 , respectively. If we applv lemma 3 to 
(27), it follows from (9) and (26) that 



^( 1 ~Wa) w ^ 1) ^( 1 +Ja) w - 



(29) 



Next we apply the corollary of lemma 1 to (28) 
and use the inequality (23). This shows that 

!Ff»|<|- (30) 

The inequalities (29) imply 

| M tf/_ F «>|<l (31) 

In fact, we have, for €<4A, 

1-e/SA- ^\A 
and, therefore, from the first inequality in (29), 

MW /_yu)<_<Lya>. (32) 

By the definition of F 8 (1) and inequality (30), 

|F 5 (1) |<2A 

If this is combined with inequality (32), one of the 
two inequalities equivalent to (31) is at hand. The 
second inequality is proved similarly. 

Finally, adding inequality (22) to (31) the desired 
formula 



\iiW-V\<e, for [ M |< 



Mo 



is proved. 

We summarize this result as a theorem. 

Theorem 2. If assumptions (A), (B), (L), (E), and 
(S) are satisfied, then the mean duration W (P) of 
the random walk satisfies the asymptotic relation 

lim M W r (P) = F(P), 

uniformly in G, where V(P) is the solution of the 
elliptic boundary value problem (11). 

V. Bounds for the Mean Duration in a 
Simple Case 

From the viewpoint of numerical application the 
preceding theory is of interest mostly as a basis for 
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finding reasonable bounds for the expected duration 

There exists several methods for this. But since 
this investigation was originally motivated by ques- 
tions connected with the "Monte Carlo' ' technique 
of solving differential equations experimentally by 
means of random walks, we shall disregard methods 
which amount themselves essentially to an approxi- 
mate solution of such a differential equation. A 
general procedure for obtaining bounds for W (P) is 
given in [2]. Here we shall limit ourselves to the 
special case that L[u] is a constant multiple of 
Laplace's operator. Then it is possible to obtain 
much more precise information, even if the space is 
again assumed to be n-dimensional. 
In other words, we require that 

PiK(P) = k.8 iK ,a i (P) = 0, i,K=l,2, • • • , ft(33) 

The differential problem (11.) now assumes the form 

AF+|=0m6?,F=0on(7 

where A is Laplace's operator in n dimensions. 
The constant k is related to the mean square of 

n 

the step length, s 2 =^ ba{P) by the formula 

i=l 

s 2 —rtixk-\-o (/x) 

There is therefore a mi>0 such that for |m|<Mi it is 
possible to use s 2 itself as a parameter /x. Without 
loss of generality we can therefore set 



--S 2 



and, hence, 



k=- 



n 



and 



AV+2n=0 in G, V=0 on C 



(34) 



(35) 



(36) 



For those arguments below, where we let n ap- 
proach infinity we need the further assumption that 
Hi can be chosen independently of n. Kandom walks 
of the type described above will be called symmetric. 

The two most interesting special cases of the type 
considered are 

(a) The random walk takes place in an orthogonal 
net of mesh length h with transition probability 1/2 n 
in each of the 2n possible directions. 

Here a i = ) biK = h 2 b iK ln. 

Assumptions (B), (L), and (E) are satisfied, if /* is 
defined by 

ix = h 2 . (37) 

Then (35) and (34) are satisfied. 

(b) The transition probability is normal with 





di=Q, 


bi K =<r 2 SiK 






(38) 


If we define 














a 


i = n<r 2 






(39) 


formulas (34) 


and (35) 


remain again va 


id. 




Theorem 3. If the 
the mean duration of 
mensional sphere of rad 
formula 


conditions (33) 
a random walk 
ius a satisfies the 


are satisfied 

in an «-<li- 

asymptotic 




W{P)=- 


a2 r 2 [i+. 


»(«)] 




(40) 



where r is the distance of P from the center, s 2 the 
mean square of the step length, and lim e(s) = 0. 

n—>oo 

The most interesting aspect of this formula is its 
asymptotic independence of the dimension. 

Proof: If G is a sphere in n dimensions, the 
solution of problem (36) is V=a 2 — r 2 . Theorem 3 
is therefore an immediate consequence of theorem 
2 and formula (34) . 

Corollary 1: Let v be the volume of the n-di- 
mensional sphere of radius a. Then the value of 
W(P) at the center of the sphere is 

W(0)=T^ n (^^^ir-' 1 s- 2 v 2/n [l+€(8)]. 

Corollary 2: 

lim W(0)/n= A(s), where A(s) ^ (41) 



This follows by a simple application of Sterling's 
formula to the right side of the preceding formula. 

If G is not a sphere, the equality (40) can be re- 
placed by an inequality, if we make use of the 
following lemma. 

Lemma 4- Let P be a fixed point in n-dimen- 
sional space and denote by {B} the set of all n-di- 
mensional domains of constant volume v containing 
P in its interior and having a boundary which is 
sectionally in Class C (2) . Let u B {P) be the value at 
P of the solution of the differential problem 

Au+1 = in #,^=0 on C. 

Then u B (P) assumes its maximum with respect to 
the Class {B}, if B is the sphere if with center at P. 
Prooj: Let r be the distance PQ and denote by 
G B (Q) Green's function in B, if one argument point 
is fixed at P. For the sphere K of radius a about 
P Green's function is of the form 



G K (Q)=c 1 (n)[r 2 - n -a 2 - n ] } n>2 



(42) 



where Ci(ri) is a constant depending on the dimension 
n. For any domain B we have 
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«s(P)= f G B (Q)dQ (43) 

Denote by v B (P), P>0 the volume of the "level" 
surface G B (Q)=p- For the sphere K we have, 
from (42), 



Mp)=c 2 (w)G>+a 2 -T'«-»> 

We now make use of the. inequality 

Vb(p)<v k (p), for all B of {B} 



(44) 



(45) 



which is an easy consequence of the classical ex- 
tremal property of the capacity of a spherical shell, 
(cf. Szego [6], in particular p. 588. Szego's proof 
for the three-dimensional case extends readily to 
higher dimensions.) From (43) we have 

u B (P) = - j P dv B (p) = — \ pv B (p) L J v B (p)d P . 

The first term in the right member is zero, for 

v B (0) = v and lim pv B (p)=0, since pv K (p) — which 

p — >00 

by (45) at least equals pv B (p)— tends to zero, as 
p->°°, in consequence of (44). A reference to (45) 
now completes the proof of the lemma for n>2. 
The modifications for n=2 are trivial. 

From this lemma and corollary 1 of theorem 3 we 
obtain immediately the following result. 

Theorem 4- The mean duration of a symmetric 
random walk in an ^-dimensional domain of volume 
v satisfies the inequality 

W(P)<TT- l T 2 ' n 0^) s- n v 2/n [l+e{s)}. (46) 

Corollary 2 to theorem 3 implies that the right 
member of inequality (46) tends to infinity with 
order n, as n-> oo f provided the volume is kept 
constant. But an explicit calculation shows that 
the increase, for small integers n is slower. In fact, 
if we write the right member of inequality (46) in 
the form 



c(n> 2/w s- 2 [l + €($)], 



we have 



n 


1 


2 


3 


4 


5 


6 


7 


c{n) 


. 250 


.318 


.385 


.450 


.515 


.578 


.642 



VI, Random Walks in an jn-dimensional 
Cube 

The n-dimensional cube is another domain for 
which problem (36) can be easily solved. It seems 
intuitively plausible that the value of W(P) at the 



center of such a cube will have the same order of 
magnitude in n as the value for the sphere, given by 
(41). This is, however, not the case. We shall, in 
fact, prove the following theorem. 

Theorem 5: For a symmetric random walk prob- 
lem let W a n denote the mean duration of walks start- 
ing at the center of an ^-dimensional sphere of volume 
v, and denote by W c n similarly the mean duration 
for a walk starting at the center of an ^-dimensional 
cube of the same volume v. Then 



but 



lim W c JW s n = 



lim n e W c n /W s n = oo ; (e>0, arbitrary). 



Proof: It suffices to consider the unit cube B, de- 
fined by \Xi\<^—, (i = l, • • -,4 The eigenfunctions 
of Au in B are 

n 

u K K (x h • • • , x n )= ncos KiirXi^KiyO.Ki odd, 

1 n 2=1 

and the corresponding eigenvalues are 



1 2=1 

These eigenfunctions are orthogonal. Their norm 
is 2 -w , since 

• • • n cos 2 KiirXi dxi • • • dx n = 2 n > 

-\ J-i*'=i 

It follows that the multiple Fourier coefficients 
Ck^ • • • K n of the representation 

1 = XUx ...K n c Kl ... K n u Kl . . . K n {x x . . . x n ), (odd K^^Kn) 
valid in B, are 

If we now substitute for u in Ai£+1 — a Fourier 
series 



u- 



K \X\, . . . , X n ) f 



(odd K u ..., K n ) 
comparison of coefficients shows that 

that is, the value of u at the center x t = of the 
cube is 

\T/ n i=l I i=l i=l 

(od&K u ...,K n ). 
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This n-fold sum can be simplified by means of the 
substitution 

co -x 2 K\ 



/ n \-l f co -x Z Kt 

(s«) -X « '■■ d 



followed by an interchange of summation and inte- 
gration, which is not difficult to justify (see, e. g. [1] 
I, p. 331). Then we obtain 



(K l9 . ..,K H odd), 
that is 



n (_1)(^,~l)/2 



K t 



e-xK?dx, 



/A\n /»oo / co r (2j+l) 2 i)» 

« (o »^- J Q J. fe^'w! *• <"> 

In order to investigate the asymptotic dependence 
of u(0) on n we make use of the formula 



y^ g-j(7«+r) : 



% X 72= —a 



exp{~xT 2 + {xT+Tin) 2 /x] 



(see, e. g. [5], p. 32), which is essentially the reciproc- 
ity formula for the ^-function. If we make the sub- 
stitution T=iy/x and integrate with respect to 
y between and 7r/4 we obtain from this formula 
the series 



/(*)=■-§+ 



Vi?r< 



exp [—(y + im) 2 /x]dy 



for the function 



Setting t= (y + irn/ -yj x in the integrals above, this be- 
comes 



/(x)= _| + v;{J o 



T/i^Jx 



e-t 2 dt+ 



(»+i)/V* 



n = lJir(»-i)/Var 



From (48) we conclude 

/(*)<-§ + V^V^^J for z>0 

and from the definition of j(x), 

j{x)<e-\ 
Hence, 

-.(o)-.-(|)-{j;+X'+fyw^| 

(!„,)■ (1 _ <)+ (i C! )'j, 



(48) 



<- 



« + 



where e is an arbitrarily small positive number and 
Ci and c 2 are constants less than tt/4. Therefore 

Um lt ( ) =0i 

7l-*co x ' 

The first part of theorem 5 is now an imme- 
diate consequence of theorem 2 and formulas (36) 
and (41). 

In order to estimate the rate of convergence of 
16(0) from below, we observe that 

u(0)> ^" 2 (^Y P j n {x)dx. (49) 

In the right member we use the inequality 

j(x)>~+J^r e-* 2 dt-^r e~ i2 dt, 

* JO JtcIWx 

which follows from (48), combining it with the 
inequality 

§° e- t2 dt<±e-**,forz>0, 

which is readily proved by an integration by parts. 
Then we obtain 



/(*)>!( 



1- 



8Vz 



rr-3/2 



x )>l(l-2x N )>~e- 2xN 



for 0<x<l; iV>0, arbitrary. Inserting this in 
(49), we have, finally, 

^(0)>7r" 2 Texp [-2nx k ]dx=Y^ m r e'^r^-'dr 

> const. n 1/N . 

This, in combination with formulas (36), (41), and 
theorem 2, proves the second part of theorem 5. 

VII. Random Walks in n-Dimensional 
Ellipsoids 

In this section we collect some results concerning 
the mean duration of random walks in 7i-dimensional 
ellipsoids. 

The "oblongness" of such an ellipsoid can be 
defined as the ratio of the longest to the shortest axis. 
The mean duration of a random walk of the class 
considered in the last two sections will be less for a 
walk starting from the center of such an ellipsoid 
than it would be for a sphere of same volume. The 
ratio between these two quantities will give an idea 
of the deterioration of the estimate (46} for oblong 
convex domains. For certain domains it may even 
be advantageous to estimate the mean duration of 
the random walk by its value in a circumscribed 
ellipsoid rather than by formula (46). We shall 
make use of a lemma concerning the arithmetic and 
geometric means of n numbers which may be of 
some interest in itself. 



933527—51- 
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Theorem 6. The mean duration of a symmetric 
random walk in an n-dimensional ellipsoid with 
semiaxes a h . . ., a n is 

W E (P) = s- 2 [l + e(s)}n(l-±§)l±± (50) 

Proof. It is readily verified that the solution of 
problem (36) in such an ellipsoid is 

/ n 1 



V(P) = n(l-±f\l±± 



This fact and theorem 2 prove formula (50). 

Theorem 7. If W S (P) is the mean duration of a 
symmetric random walk for a sphere having the 
same volume as the ellipse, then 



[!+€(«)] {n ar 2 



W E (0)/W s (P) 
Proof: The volume of the ellipsoid is 



i=i 



(51) 



(52) 



If we divide formula (50) for P=0 by the formula 
in corollary 1 to therorem 3, and eliminate v by 
means of formula (52), we obtain formula (51). 

From now on we shall assume that the oblongness 
of the ellipsoid is prescribed. In view of the remarks 
at the beginning of this section it is of interest to 
investigate the minimum, of the right member in 
(51) under this restriction. To this end we prove 
the following lemma: 

Lemma 5: Let tj (j=l, . ..., ri) be n positive 
numbers such that the ratio of the largest to the 
smallest has a prescribed value p. Then the 
quotient of the geometric and the arithmetic mean 
of these numbers will be smallest when 

H == *2 =r ' ' * =z t / r) *r+l =z fr+2 ==: * " * == ^n.) 

where r is an integer that satisfies the inequalities 

x(p)n<r< x (p)n+l, 



with 



x(p)=- 



" p — 1 log p 
Proof: Without loss of generality we can assume 



We have to minimize the quantity 



Fth, 



., tn-l)~- 



-- n «,/W 



(53) 



(54) 



under the side conditions (53). As a function of t s 
alone, (l<s<7i), F has its only stationary value 
when 



y^ f tj—nt s = 



;'=i 



i.e., 



t s =- 



1 






This stationary value is a maximum, since F is zero 
when t s = or t s = co , Hence, as a function of t s 
alone, F has no relative minimum and it assumes 
therefore its minimum for given values of tj(j-{-s) 
either when t s =t 8 - 1 or when t s =t s+1 . It follows that 
at the minimum of F, the tj consist of two groups 
such that 

ti=t 2 = • • • =t T ,t r+l = • - • =t n ^ 1 = t n , \<r<n—l 

the minimum will therefore be among the values of 

p n-r 



*(r) = 



[r+(n-r)p\* 



1 <r<n— 1. 



We find that \f/(r) has exactly one stationary value, 
which occurs at 



r=r =n x (p)=n 



p - 1 log 



g P/ 



This stationary value is a minimum, for \f/(r) is con- 
tinuously differentiable and positive in the interval 
— oo <><^p/(p — \) it is infinite at r = — oo and r= 
np/(p— -1), and r Q lies in this interval. Moreover, 

lim x(p) = o> x(p) increases with p, and x(°°)=l. 

p-»i ^ 

Hence 

n 

- <r <>, 



If r is an integer, this is the admissible value of r 
for which \f/(r) assumes its minimum. Otherwise, it 
is one of the two neighboring integers. This com- 
pletes the proof of the lemma. 

Corollary 1: For sufficiently small p>l the in- 
teger r of lemma 5 differs by less than unity from 
n/2. As p increases, so does r. For sufficiently 
large p we have r=n— 1. Also, r is never less than 
n/2. 

Corollary 2: For n=?> we have always r—2. 
Proof: One shows easily that, for p>l, 

*(i;#(2)>i. 

Theorem 8: With the notations and assumption 
of theorem 7 denote by m(co) the minimum of the 
ratio W E (0)/W s (0) for all ellipsoids of given oblong- 
ness CO. 
Then 



lim {m(w).[l +€(«)]} = 

[x(V«)w 2(1 -* ( V5» + co 2 (l - X ( Vw))] - 1 . 



(55) 
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Prooj: Without loss of generality we set 

d\>(^2> • • • >dn, 



hence 






We apply lemma 5 to formula (51) with 



tj=a- 



1. e., 



-=(if)=' 2 - 

The minimum of the right side of (51) is then 
[l + e(a)]na^ 2r/n a- 2in - r)/n /ra^ 2 + (n-r)a- 2 == 

[!+€(*)]. (56) 



— CO +1 

n n 



As n-> oo, r/n->x(p) = x(iju ) ), and formula (55) is 
at hand. Theorem 8 shows that for sufficiently 
large n and co> 1 



m(co)<^co' 



2(X (V ^)-1) 



<1 



On the other hand, m(co) does not tend to zero as 
n— > oo . For large n, random walks starting from the 



center of any ellipsoid of prescribed oblongness are 
on the average longer than those in n-dimensional 
cube of same volume. If r is kept constant for all n 
in (56), it represents the ratio W e (0)/W s (0) for 
ellipsoids of given oblongness r of whose axes have 
maximal length while the remaining n—r have 
minimal length. For this type of ellipsoids 
W E (0)/W s (0) tends to 1 as n-^> co . The correspond- 
ing random walks are, therefore, on the average 
longer, for large n, than for ellipsoids of the minimal 
type. 
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